Critical-current density from magnetization loops of finite high-T c superconductors 
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We analyze the effects of demagnetizing fields in the magnetization hysteresis loops of type-II super- 
conductors, by a model that allows the calculation of current and field profiles in finite cylindrically 
symmetric superconductors in the critical state. We show how the maximum in the magnetization 
curve gradually shifts from negative applied field values to the central position with decreasing sam- 
ple thickness. From the analysis of the calculated field profiles, we demonstrate that one can obtain 
the intrinsic field-dependent critical-current density of the superconductor with great accuracy by 
measuring the magnetic response of superconductors with large aspect ratio with the applied field 
parallel to the shortest dimension. 
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The critical current density J c is one of the key param- 
eters of high-Tc superconductors. In particular, know- 
ing the dependence of J c on the internal magnetic field 
Hi = B I no in the superconductor is a key concern for 
the study of the current-carrying states of the supercon- 
ductors. The J c (\Hi\) dependence is obtained from two 
general types of experiments: electrical and magnetic. In 
electrical transport measurements, the reliability and re- 
peatability of the measured values of the critical current 
is reduced because of the difficulties in making contacts in 
the sample and in choosing a voltage criteria for defining 
J c . Moreover, in general large currents circulate through 
the sample during the measurements, so that they pro- 
duce a magnetic field that modifies the applied field in 
a inhomogeneous way. This makes the estimation of the 
intrinsic J c (\Hi\) dependence difficult. The other general 
method for obtaining J c is from magnetization measure- 
ments, typically isothermal magnetization loops M(H a ), 
where H & is the applied field. J c is obtained from the 
width AM of the M(H a ) hysteresis loop using 



J C (H) = AM(H)/d, 



(1) 



where d is a length characteristic of the sample size and 
geometry; for a cylinder of radius R, 



J C {H) = 3AM(H)/2R. 



(2) 



This relation follows from the critical-state model jjj. 
The magnetic method also involves some approximations. 
First, the method was originally derived for the case of 
J c independent of \Hi\ (I). If instead J c depends on \H{\, 
as in most actual superconductors, what is obtained from 
the width of the loop is not the required intrinsic J c (|-ffi|) 
function, but a different function, J c (\H a \). These func- 
tions are approximately equal only if some extra condi- 
tions concerning internal field homogeneity are met, as 
discussed in and also below. The second important 
approximation is that the method is only theoretically 
justified for infinitely long samples. The validity of the 
method for extracting J c from the width of the loop is 
therefore questionable in the realistic case of finite super- 
conductors. 



Calculations of the magnetic response of finite super- 
conductors in the critical state including demagnetizing 
effects have been recently presented for strips || and 
cylinders j7|-|j|, following previous works on very thin 
strips jlo| and disks jllj. Very recently, Shantsev et al 
]l2t have discussed important features of the effects of de- 
magnetizing fields in the hysteresis loop. In particular, 
they demonstrated theoretically and experimentally that 
for very thin strips in perpendicular field the peak that 
appears in the reverse magnetization curve is located not 
at negative fields, but at the central position. 

In this work, we systematically study the effect of de- 
magnetizing fields in the magnetic response of finite su- 
perconducting cylinders in the presence of a uniform ap- 
plied field, and discuss the implications in J c extraction. 
We first introduce a model which allows us to compute 
current and field profiles and magnetization loops of su- 
perconducting cylinders with the same intrinsic parame- 
ters but different aspect ratios, in order to study in detail 
the effects of demagnetizing fields. These will be analyzed 
in relation to two particular features: the position of the 
peak in the magnetization curves as a function of the 
sample aspect ratio, extending the work of Shantsev et 
al |l2j , and the relation between the shape of the hystere- 
sis loop and the field dependence of the critical current 
density for the different aspect ratios. 

Our model, which simulates the process of penetration 
of supercurrents inside a superconductor, is based in the 
fact that any modification of the applied field results in a 
change in the superconductor current distribution in or- 
der to minimize the change in the magnetic energy. The 
current distribution in the initial magnetization curve 
(after the sample has been zero-field cooled and a mag- 
netic field is applied) can therefore be obtained from the 
magnetic energy minimization. We will obtain the re- 
verse current distributions (once the field has reached a 
maximum value and is decreased) by the conventional 
procedure of superposing a current distribution with op- 
posite sign to the 'frozen' field profiles as in |l0|. In all 
cases, we assume that B — /io-ff, which holds for fields 
H c i « H « H C 2, as is usually done in critical-state 
modeling. 
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Consider a cylindrical type-II superconductor of radius 
R and length L located in a uniform applied field, H a , 
directed along its axis. We use common cylindrical co- 
ordinates (p,9,z), z being the direction of the axis of 
the superconductor. Owing to the symmetry of the sys- 
tem (which reduces the problem to a two-dimensional 
one), all supercurrents flow with angular direction. We 
divide the superconductor in a regular grid of n x m 
coaxial rings in which linear currents can flow. The 
magnetic flux that threads one of these linear circuits 
at the position (p, z) due to the external applied field 
is z) = poH a irp 2 , while the flux that threads the 

same circuit due to all currents circulating in the su- 
perconductor is <& l (p, z) — J2 P > z ' M (Pt z , p', z ')I(p', z '), 
where I(p', z') indicates the current in the circuit (//, z'), 
and M(p, z, p' , z') the mutual inductance between the cir- 
cuits (p' , z') and (p, z). The self-inductances M(p, z, p, z) 
are calculated from the mutual inductance between two 
close linear circuits ||. The model details are described 
in Refs. |^,|l3|]. In this paper we have used typical values 
of n x m ~ 120 x 20 for thin samples and n x m ~ 60 x 60 
for samples with larger L/R ratios. 

Let us assume we have a given current distribution 
corresponding to an applied field H a (if we are calcu- 
lating the first point after the initial state, then the 
initial current distribution is zero everywhere). Set- 
ting a current I at a circuit (p, z) requires an energy 
E = I$ l (p, z) while it contributes to reduce the energy 
(current has opposite sign to H a ) by a quantity IQ a (p, z). 
We find in this way the circuit that yields the largest de- 
crease of energy and set a current / there. The new 
currents are set accomplishing the chosen material law 
J = J c (\Hi\), \Hi\ being the modulus of the total field 
Hj, and / = J c (\Hi\)(RL/nm). When no new currents 
minimize further the energy, we calculate the magnetic 
field inside the superconductor and change the value (not 
the distribution) of the already induced currents to ac- 
complish the material law. At this point, we calculate the 
magnetic moment resulting from all the circulating cur- 
rents and all the other relevant magnitudes. After this, 
the applied field can be increased again, and the pro- 
cess is restarted from the existing current distribution. 
The reverse stage (corresponding to decreasing H a from 
ff max to — H max ) is calculated by superposing the frozen 
current penetration set at H max to the one induced in 
the reverse stage (calculated in a similar manner as de- 
scribed for the initial curve). This procedure is typical 
of the critical state model [n0[ . 

Our model allows the implementation of an arbitrary 
J c (\Hi\) dependence. The following discussions and con- 
clusions are valid for any dependence as long as J c is a de- 
creasing function on which is physically reasonable. 
For illustrating our results, we choose an exponential de- 
pendence J c = J c o exp(—\Hi\/Ho), where J c o and Hq are 
positive constants. The exponential dependence has been 
successfully applied to high-T c superconductors [ ]l4| , |l5| |. 
A useful parameter for the analysis is p — J c qR/Hq, 



which takes on large values for a strong dependence of 
J c on | Hi | and tends to for J c independent of field [fli) . 
Realistic values of p for high-T c superconductors range 
from 1 to 10 Jy|. 

In Fig. 1 we show the calculated M(H a ) curves for 
different values of L/R, for the cases p = (Bean's 
model), 3, and 10. M and H a are normalized to H p = 
-ffoln(l +p), which corresponds to the penetration field 
for an infinite cylinder. We find that the calculated 
M(H a ) loops for sufficiently large samples coincide be- 
tween our numerical accuracy with the known analytical 
results for infinite cylinders jjyj. Let us now discuss the 
differences observed in the loops for each value of p, which 
are only due to the sample geometry. In all cases, the de- 
magnetizing field enhances the initial slope of both the 
initial and reverse curves. We have analyzed this effect 
in detail, observing good agreement with experimental 
data measured for niobium cylinders of different lengths 
in 1 13 1; we obtain the same values for the initial slope of 
the M(H a ) curves as those calculated by Brandt and 
Chen et al Q for a wide range of L/R values, with less 
than 1% deviation. The shape of the central part of the 
loop is not changed in the Bean case (p = 0), since even 
in a finite sample the width of the loop is proportional to 
the the constant value of J c . However, important vari- 
ations are observed for non-constant J c . When p ^ 0, 
a peak appears in the reverse magnetization curve. In 
Fig. 2 we show the calculated dependence of the peak 
position -ffpcak on L/R for p=3 and 10. 

The general trend of the peak position is that it tends 
towards H a = with decreasing sample L/R ratio. A 
similar tendency was found in |17| for a three parame- 
ter J c (Hi) function and using an iterative method valid 
only for the fully penetrated stage. Our calculated re- 
sults give the correct known limit when the sample is 
very large (solid line in Fig. 2 indicates the limit for 
infinite cylinders 0|). Our results are also compatible 
with where it is said that the peak tends towards 
H a = for thin samples without reaching exactly the 
central position (this is in contrast with the strip case, 
for which the peak position was shown to be zero for very 
thin strips @,||). 

The peak position and the other observed features of 
the hysteresis loop can be understood as follows. We 
start by discussing the situation in the simpler case of 
an infinitely long superconductor with a given J c (|_Hi|) 
function. After reaching the maximum applied field 
H max , when we reverse the magnetic field sweeping di- 
rection in the hysteresis loop, supercurrents are induced 
at the surface of the cylinder with a direction opposite to 
the shielding currents of the initial magnetization curve. 
These reverse supercurrents gradually enter the sample, 
confining the original currents to the interior. When H a 
decreases, the internal field in the region penetrated by 
reverse currents also decreases, these currents become 
larger, and therefore the magnetization increases. Pass- 
ing through H a = 0, the internal field becomes negative 
at the surface regions of the superconductor. Further de- 
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creasing H a makes the Hi profile negative, so that \Hi\ 
again increases. Since J c decreases with the absolute 
value of Hi, the currents become low and M decreases. 
Then, a peak appears in the magnetization at some neg- 
ative value of H a , for which an averaged value of |ffi| is 
minimum. This process is depicted in Fig. 3a, where we 
show the calculated field profiles in the cylinder midplane 
corresponding to the case p = 10 for a long cylinder with 
L/R = 10, which schematically represents the infinite 
case. 

There is a key difference for this case and the behav- 
ior observed for a realistic finite sample. Whereas in an 
infinite sample the internal field at a given point in the 
superconductor has only contributions from H a and the 
field created by the currents exterior to this point HQ , 
there is contribution from all the circulating currents at 
all points in a finite superconductor. This self-field con- 
tribution produces an effect in the central region of the 
loop (\H a \ small) similar to shifting the applied field up- 
wards to higher values for H a > and downwards to 
lower values for H a < 0. This explains the increase in 
the initial slopes of the loops shown in Fig. 1, as well as 
the appearance of a peak in M for applied field values 
H a much closer to the central position. We show in Fig. 
3b the calculated field profiles for the same case, p = 10, 
as in Fig. 3a, but for a thin disk with L/R = 0.1. The 
large effect of demagnetizing fields is made manifest by 
the almost constant field profiles characteristic of high 
applied fields already achieved for low values of | H a \ . 

This has important consequences for the method used 
to extract J c (\Hi\) from the width of the magnetization 
loop M(H & ). As discussed by Chen and Goldfarb ||] 
there are two requirements for using the method: (1) The 
magnetization on ascending and descending branches of 
the hysteresis loop at a given field H a must correspond to 
fully penetrated states; and (2) the maximum deviation 
of J c (\Hi\) in the sample from the value of J c (\Hi\) for 
Hi = H a must be small. Condition (1) is easily fulfilled 
as long as one takes care to measure a reverse magneti- 
zation curve starting from a sufficiently large maximum 
applied field value. The key condition is then the sec- 
ond, which is understandable taking into account the 
fact that the formula for J c extraction is based on the 
Bean's model for constant J c so large inhomogenieties in 
the local currents will yield wrong results. Comparing 
our results for field profiles calculated for the same ma- 
terial but with different dimensions, see Figs. 3a and 3b, 
one can clearly see that the field profiles in the case of 
thin samples are much more spatially uniform than in 
the case of long superconductors. (Actually, the value 
of the internal field for thin samples is basically equal 
to the external H a value except within a narrow region 
around |_ff a | = 0.) There can be several reasons for that. 
Currents are low, except when |£f a | is very small, ow- 
ing to the contribution of the demagnetizing fields to Hi 
in the J c (\Hi\) function, producing small contribution to 
the total field. Moreover, in a thin disk a much smaller 
total current is flowing as compared with a bulk cylin- 



der. Therefore, thin samples in transverse geometry are 
an optimum case for obtaining J c from the width of the 
magnetization loop. This conclusion is confirmed in Fig. 
4, where we compare the J c (H a ) function extracted from 
the width of the loop with the intrinsic J c (\Hi\), for the 
cases p = 3 and 10, and for a long sample {L/R = 10) 
and a thin one (L/R = 0.1) for each case. The agree- 
ment between both J c functions is clearly better for the 
thin samples than for the long ones for all magnetic field 
values, and particularly for the low-field region which in- 
volves in the long sample a large measurement error. Our 
model also shows that for higher values of p the agree- 
ment between J c functions is better, which implies that 
measuring thin superconductors is specially important 
when the expected dependence of J c on \Hi\ is strong. 

We conclude from this analysis of demagnetizing effects 
in realistic finite superconductors that the extraction of 
the critical current density from magnetization measure- 
ments is best done in thin film geometry in perpendicular 
field, or, if thin films are not available, with the field ap- 
plied along the sample shortest dimension. In this way, 
the superconductor critical current density and its depen- 
dence on the internal field can be precisely obtained by 
measuring the hysteresis loops of superconductors using 
the non-destructive, reproducible, and widespread mag- 
netization measurements. 
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FIG. 1. Magnetization loops for the exponential depen- 
dence of the critical current, for (a) p = (Bean model), 
(b) p = 3, and (c) p — 10, and for different values of the 
length-to-radius, L/R = oo (solid), 1 (dashed), and 0.1 (dot- 
ted). 

FIG. 2. Calculated applied field at which a peak in the 
reversal magnetization of superconducting cylinders occurs, 
flpcak, as a function of L/R. Filled circles correspond to p = 3 
and open circles to p = 10. Horizontal lines represent the 
known values for infinite cylinders. 

FIG. 3. Internal magnetic field Hi profiles in the midplane 
of the superconductor cylinder for the case p = 10 and L/R = 
(a) 10 and (b) 0.1, corresponding to the reverse curve after a 
maximum applied field of H max — 1.5H P . The values of the 
applied field ff a range from 1AH P to -1.4H P in steps of 0.2H P 
(from top to bottom). 

FIG. 4. Critical current density calculated from the width 
of the magnetization loop AM using Eq. (2) for L/R— 0.1 
(dashed line), and for L/R — 10 (dotted line), and from 
the analytical expression of the exponential dependence (solid 
line). The case (a) corresponds to p = 10 and (b) to p = 3. 
H a is the external applied field at which AM is evaluated, 
whereas H\ is the internal field upon which J c depends. 



4 



0.5 n — ' — i — ' — i — ' — i — ' — i — ' — i — ' — r 



k p=o 






k (a) ^ 



-0.5 H — i — I — i — I — i 





L/R 



